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ABSTRACT 

Recent studies by a number of independent collaborations, have correlated the CMB 
temperatures measured by the WMAP satellite with different galaxy surveys that 
trace the matter distribution with light from the whole range of the electromagnetic 
spectrum: radio, far-infrared, optical and X-ray surveys. The new data systematically 
finds positive correlations, indicating a rapid slow down in the growth of structure in 
the universe. Individual cross-correlation measurements are of low significance, but we 
show that combining data at different redshifts introduces important new constraints. 
Contrary to what happens at low redshifts, for a fixed f2 m , the higher the dark energy 
contend, the lower the ISW cross-correlation amplitude. At 68% confidence level, 
the data finds new independent evidence of dark energy: fl\ = 0.42 — 1.22 . It also 
confirms, to higher significance, the presence of a large dark matter component: f2 OT = 
0.18 — 0.34, exceeding the density of baryonic matter, but far from the critical value. 
Combining these new constraints with the prior of a flat universe, or the prior of an 
accelerating universe provides strong new evidence for a dark cosmos. Combination 
with supernova data yields f^A = 0.71 ± 0.13, f2 m = 0.29 ± 0.04. If we also assume 
a flat universe, we find fl\ = 0.70 ± 0.05 and w — —1.02 ± 0.17 for a constant dark 
energy equation of state. 



1 INTRODUCTION 



In the last few years a new cosmological scenario with a sig- 
nificant smooth Dark Energy (DE) component has emerged. 
The Cosmic Concordance Model (CCM, from now on) is a 
spatially flat universe with baryons (fib ~ 4%), cold dark 
matter (Qcdm ~ 23%) and a significant DE component 
(Qa ~ 73%). The model is well supported by the supernova 
type la observations (SNIA) (Riess et al. 1998; Perlmut- 
ter et al. 1999), observations on large scale structure (LSS) 
(Tegmark et al. 2004; Percival et al. 2001) and the cos- 
mic microwave background experiments (CMB), in particu- 
lar by the recent WMAP experiment (Bennett et al. 2003) . 
The energy density of the universe seems dominated by the 
unknown DE component, presenting a formidable observa- 
tional and theoretical challenge. The three key observational 
probes measure complementary aspects of the cosmological 
parameter space. The SNIA indicate that the universe is 
accelerating but present data is degenerate for alternative 
cosmological scenarios. The LSS observations constrain Q m 
but leave the DE question unanswered. Constraints from pri- 
mary anisotropics in the CMB indicate that we live in a flat 
universe but require a prior on the value of the local Hub- 
ble rate _ffo(Blanchard 2003). Assuming that the universe is 
well described by a ACDM model, combining all these three 
observations gives us the cosmological CCM model. 



The Integrated Sachs Wolfe effect, ISW, (Sachs & Wolfe 
1967) is a direct probe for the (linear) rate of structure for- 
mation in the universe. Secondary anisotropies in the CMB 
appear because of the net gravitational redshifts affecting 
CMB photons that travel through an evolving gravitational 
potential $. These secondary temperature anisotropies are 
therefore correlated with local, evolving, structures on large 
scales. The correlation is negative when structures grow, as 
increasing potential leaves a cold spot in the CMB sky, and 
positive otherwise. In a flat universe without DE (Einstein- 
deSitter, or EdS, model) this cross-correlation is expected 
to be zero because the gravitational potential remains con- 
stant, despite the linear growth of the matter fluctuations. 

The rate of structure formation in the universe can also 
be measured by galaxy peculiar velocities or galaxy redshift 
distortions, on very large scales through the so-call /3 param- 
eter determination (Peacock et al. 2001; Pope et al. 2004). 
The ISW effect provides an independent and complemen- 
tary probe of the same effect. Independent, because it uses 
temperature anisotropies instead of the velocity field, and 
complementary, because of the different assumptions and 
systematics that relate measurements with theory. Despite 
recent advances in the size of galaxy redshift surveys such 
as SDSS and 2dFGRS, the spectrum of matter fluctuations 
P(k) oc< 5(k) 2 > is quite difficult to measure directly over 
very large scales (Tegmark et al. 2004; Percival et al. 2001; 
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Gaztanaga & Baugh 1998). Part of the problem is that 
matter correlations fall quickly to zero on scales larger than 
30 Mpc/h (k < 0.1 h/Mpc). In contrast, fluctuations in the 
gravitational potential go as $(fc) cx S(k)/k 2 and therefore 
extend over larger distances, which makes the signal more 
detectable (see also comments to Fig. QJ. The ISW cross- 
correlation traces the gravitational potential, "!>, and thus 
provides a new window to study the largest structures, ex- 
tending over several degrees in the sky or tens of Mpc/h at 
the survey depth. 



2 GROWTH OF DENSITY PERTURBATIONS 

Gravitational evolution of matter fluctuations, 5 = p/p — 1, 
is dependent on the cosmological model via the evolution of 
the scale factor a — a(t). Compared to a static background, 
a rapidly expanding background will slow down the collapse 
of an over dense region. In the linear regime, a small initial 
perturbation <5o grows according to the growth factor D(t): 



6{t) = D(t) S 



(1) 



which, under quite generic assumptions, eg (Gaztanaga . 
Lobo 2001; Multamaki, Gaztanaga & Manera 2003; Lue . 
Starkman 2004), follows a simple harmonic equation: 



d 2 D 
dr/ 2 



+ 2 



H 

Jr 2 



dD 



+ 3ci D = 0, 



(2) 



where rj = ln(a) is the conformal time and H = H(p) = a/a 
is the background Hubble rate (d and H are proper time 
derivatives). For a flat cosmological model with a generic 
dark energy equation of state 



p = w(z)p 
we then have: 



H 2 



fl m (l + z 



(3) 



(4) 



where Oa and Q m are the dark energy and dark mat- 
ter densities today in units of the critical density p c 
:> H 2 /(8nG). And ci is given by: 



Cl = -2 



(•>) 



H 2 (z) 

In this paper we study two cases. A generic (not neces- 
sarily flat) ACDM model where DE density is constant over 
the evolution of the universe (w — —1); and a flat ACDM 
model with a constant equation of state parameter. For those 
models we have 



ci = -(1/2) n m /(n m + n A a 3 ) 



(6) 



One may choose to compare the results to the EdS 
model: Qa = 0, fl m = 1, in which case the solution to 
Eq. © is D tx a. This means that 8 grows linearly with 
the scale factor, S tx a, while the corresponding gravita- 
tional potential fluctuation, $ ~ S/a, remains constant as 
the universe expands. For non EdS models <E> would change 
during the expansion of the universe which would turn into 
a galaxy-CMB temperature cross-correlation signal. 



2.1 The ISW effect 

ISW temperature anisotropics are given by (Sachs & Wolfe 
1967): 



A isw {fi) = T(n) - T 
To 



-2 / dz ^-(n, z) 

dz 



(7) 



where $ is the Newtonian gravitational potential at redshift 
z. One way to detect the ISW effect is to cross-correlate tem- 
perature fluctuations with galaxy density fluctuations pro- 
jected in the sky (Crittenden & Turok 1996). On large linear 
scales and small angular separations, the cross-correlation 
«4g W (#) =< A^ SH/ (ni)(5 G (n2) > is (Fosalba & Gaztanaga 
2004): 



ISW /Q \ 

w TG (8) 



g{k6) 

Wisw(z) 
Wa(z) 



2tt 



dk 



P(k)g(k9) 



dzW IS w{z)W G (z 



-Jo(kr A e) 



= 3Q m (H /c) 



■ d[D(z)/a] 



dz 



(8) 



= b(z)<t> G (z)D(z), 



where Jo is the zero order Bessel function, <j>G is the sur- 
vey galaxy selection function along the line of sight z and 
ta = ta{z) the comoving transverse distance. The power 
spectrum is P(k) = A k" s T 2 (k), where n s ~ l 1 and T(k) 
is the ACDM transfer function, which we evaluate using the 
fitting formulae of Einseintein & Hu 1998. We make the com- 
mon assumption that galaxy and matter fluctuations are re- 
lated through the linear bias factor, 5a{n, z) = b(z)S m (h, z). 

For the ACDM case the ISW effect is non-zero, and 
the kernel Wisw can be well approximated by Wisw(z) = 
—3£l m (Ho/c) 2 D(z)(f — 1), where f is the relative growth 
factor, / ~ 0, rn (z) 6 ^ 11 . Wisw decreases as a function of 
increasing redshift and goes to zero both for Q m — > and 
for Sl m — » 1. At low redshifts, the ISW effect is larger for 
larger values of Sl m , but the redshift evolution depends on 
the curvature (ie how quickly the H and D evolve to the EdS 
case). This is illustrated in Fig0 which shows how Wisw 
depends on z for different values of and il m . At high 
redshifts, the lower the value of Qa (for a fixed fl m ) the 
larger the ISW amplitude. 

In Figures and we also shown for a given flat cos- 
mology model the dependence of the Wisw on redshift and 
on the equation of state parameter w. For a given redshift 
and Q m there exists a maximum of Wisw around w — —0.5. 
This maximum would translate into a maximum in the cross- 
correlation signal wtg- If data turns out to be greater than 
this maximum this would clearly disfavor models with con- 
stant equation of state. 



Throughout the paper we made the assumption of scale invari- 
ant primordial fluctuations (n s ~ 1). For other possibilities see, 
eg, (Barriga etal 2001). 
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Figure 1. Redshift dependence of Wisw(z) in Eq.|S] for different 
values of Q m and Q\. Bottom left, top right and top left panels 
shows a fixed Q m = 0.5, f2 m = 0.3 and O m = 0.1 respectively. In 
all cases: Q\ = 0.0 (dotted blue line) , f2y\ = 0.7 (continuos black 
line) and Q\ = 1.0 (dashed red line). Bottom right panel shows 
a fixed Q\ = 0.7 and H m = 0.3 (continuos black line), Q m = 0.5 
(dotted blue line) and Qa = 0.1 (dashed red line). 
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Figure 3. Dependence of Wisw on the equation of state pa- 
rameter w for fiat models at different redshifts and values of f2 m . 
f! = 0.3 (back continuous lines), f2 m = 0.25 (red dashed lines), 
n m = 0.2 (green dotted lines). Redshifts are 0.1, 0.25, 0.5, 1.0, 
1.5 from top to bottom. 




Figure 2. Redshift dependence of Wisw m eq |U for flat mod- 
els with constant equation of state. Right panel shows a fixed 
f^m = 0.2 and left panel Q m = 0.3. In both cases w = —2 (black 
continuous line), w = —1.5 (red dashed line), w = —1 (green dot- 
dashed line), w = —0.5 (blue doubledot-dashed line), and w = 
(brown dotted line) 



2.2 Bias Self-calibration 

Linear bias is used to study how well light traces the un- 
derlying statistics of linear matter fluctuations. On these 
very large scales, fluctuations S are small and linear theory 
works very well both for biasing and gravity. We remove the 
effects of biasing in our parameter estimation by compar- 
ing the observed galaxy-galaxy correlation wgg, in the very 
same samples used for the cross-correlation, to the matter- 



matter correlation w mm predicted by each model (Fosalba, 
Gaztahaga & Castander 2003). The effects of bias are also 
redshift dependent, but given a galaxy selection function 
4>g(z), picked at z = z , we approximate the bias with a 
constant b — b(z) for that particular survey. We then have: 
wtg = b(z)wTm and wgg = b 2 (z)w mm , so that an effective 
linear bias b can be estimated as the square root of the ratio 
of galaxy-galaxy and matter-matter correlation functions: 



b = 



WGG 



(9) 



Such prescription has been shown to work well in a variety 
of galaxy models (eg see (Berlind, Naratanan & Weinberg 
2001)). The values of iVmm can be computed similar to JSJ 
by 



m = 



g(k6) = 
Wm(z) = <I>g{z)D{z) 



— / dkkP{k)g{kff) 



dzWl{z)^-J {kr A i 



(10) 



where the only difference between W m and Wg is the bias 
factor b(z) in Eq.|H]. Note how the estimation of b in Eq.[5| 
depends on the normalization of the power spectrum in 
Wm m . We choose to normalize each model by fixing ag. To 
make our results independent of this normalization we will 
marginalize over as and h. Taking flat priors and ranges 
as = 0.8 - 1.0 and h = 0.72 - 0.77. We compare the pre- 
dictions with the observational data wtg normalized to the 
CCM model bias, ie totg/6, where b is estimated from Eq. 
© using Wmm in the CCM model. Consequently, for other 
models, we will need to renormalize each of the theoretical 
predictions to the CCM model bias using a "relative bias": 
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WTG d / b = b r w™ d , where b 2 = w mm /w^ is the ratio of 
the concordance model prediction to the one in the corre- 
sponding model. We choose to estimate this relative bias at 
R = 8 Mpc/h, but the actual number has little effect in our 
final conclusions. 



3 OBSERVATIONAL DATA 

Recent analysis by independent collaborations, have cross- 
correlated the CMB anisotropics measured by WMAP with 
different galaxy surveys. The median galaxy redshifts ex- 
pand over a decade (ie 0.1 < z < 1.0) and trace the matter 
distribution with light from the whole range of the electro- 
magnetic spectrum: radio, far-infrared, optical and X-ray 
surveys (see Table EJ. The cross-correlation and error esti- 
mation techniques used are also quite different but they yield 
comparable results over the scales of interest. Compare for 
example the Montecarlo errors to jackknife errors in Fig. 3 
in Fosalba & Gaztanaga (2004). In our compilation of the 
different data sets, we average the results on fixed angular 
scales around 9 = 6°. This corresponds to proper distances 
of ~ 25Mpc/h at 0.1 and ~ 100 Mpc/h at z ~ 1.0 in 
the CCM model and avoids possible contamination from the 
small scale SZ and lensing effects, eg see Fig. 3 in Fosalba, 
Gaztanaga & Castander (2003). 

Radio galaxies from NVSS (Condon etal 1998) and 
hard X-ray background observed by HEAO-1 (Boldt 1987), 
have been cross-correlated with WMAP data (Boughn & 
Crittenden 2004a; Boughn & Crittenden 2004b), to find 
a signal of 1.13 ± 0.35 times the CCM model prediction at 
z ~ 0.9. The different biases for X-rays, b 2 = 1.12, (Boughn 
& Crittenden 2003) and for radio galaxies, b = 1.3 — 1.7, 
(Boughn & Crittenden 2002) have been taken into account. 
A compatible signal has also been found with the NVSS data 
by the WMAP team (Nolta et al. 2004). 

The cross correlation of WMAP with galaxies (17 < 
bj < 20) in the APM Galaxy Survey (Maddox et al. 1990) 
(covering about 20% of the South Galactic Cap, SGC) was 
found to be wtg = 0.35 ±0.13 \iK at scales 9 — 4— 10° with 
b ~ 1 (Fosalba & Gaztanaga 2004). The cross-correlation 
of WMAP with the SDSS DR1 (Tegmark et al. 2004) (cov- 
ering about 10% of the North Galactic Cap, NGC) have 
been done for several subsamples (Fosalba, Gaztanaga & 
Castander 2003). The first sample (z ~ 0.3) contains ~ 5 
million objects classified as galaxies in SDSS (with r < 21 
and low associated error). For this sample, which has b ~ 1, 
wtg = 0.26 ± 0.13 fiK at scales = 4- 10°. The high red- 
shift sample (z ~ 0.5) has wtg = 0.53±0.21 fiK and b 2 ~ 6. 
The SDSS data has also been cross-correlated with WMAP 
by the SDSS team (Scranton et al. 2003) using nearly 25 
million galaxies in four redshift samples. Their results are 
similar with those obtained earlier by Fosalba, Gaztanaga 
& Castander (2003) but no bias from galaxy-galaxy auto 
correlation function is given. The infrared 2MASS Galaxy 
Survey (Jarret et al. 2000), with z ~ 0.1, show a WMAP 
cross-correlation of 1.53 ± 0.61 times the CCM prediction, 
with a bias of b = 1.18 (Afshordi, Loh & Strauss 2004). 

We have selected independent measurements for which 
the bias CCM b (from wgg) is known, so that we can applied 



2 


w TG /b 


b 


catalog, Band 


0.1 


0.70 ±0.32 


1.1 


2MASS, infrared (2/im) 


0.15 


0.35 ± 0.17 


1.0 


APM, optical (bj) 


0.3 


0.26 ± 0.14 


1.0 


SDSS, optical (r) 


0.5 


0.216 ± 0.096 


2.4 


SDSS high-z, optical (r+colors) 


0.9 


0.043 ± 0.015 


1-2 


NVSS+HEAO, Radio & X-rays 



Table 1. Observed cross correlation WTc/b (averaged for 9 ~ 
4 — 10°.) of WMAP anisotropics with different catalogs. Errors 
in WTG/b includes 20% uncertainty in b. Errors in the median 
redshift z are about 10% . 



the bias "self-calibration" proposed in section H2.2I The data 
is summarized in Table and displayed in Fig.l. In the re- 
sults below we also include a 10% uncertainty in the median 
redshift. We chose the values of NVSS±HEAO-l quoted by 
(Boughn & Crittenden 2004b) as representative of both 
the Nolta et al. (2004) and Boughn & Crittenden (2004a) 
analysis. For the SDSS, we chose the values in Fosalba & 
Gaztanaga (2004) where the CCM bias b is estimated using 
Eq.||§J- Note how the selected samples are complementary. 
The samples which have large sky overlap (eg 2MASS and 
NVSS±HEAO-l) have negligible redshift overlap. When the 
redshift overlap is significant (ie in 2MASS-APM or SDSS- 
NVSS could be up to 20%) the sky overlap is small (less than 
10%). Consequently, the different samples in Table have 
less than 1% volume in common. This is negligible, given 
that individual sampling errors (which are proportional to 
volume) are of the order of 30%. 

The most significant detection in Table seems to be 
the one quoted by Boughn & Crittenden (2004b) for the 
NVSS±HEAO-l samples. Given the systematic uncertain- 
ties involved in the bias and selection function of both of 
these samples, we have checked that our results do not 
changed much (less than 20% in the area of the contours in 
Fig. 3) when we double the quoted errorbar. Doubling this 
errorbar corresponds to an additional 50% systematic uncer- 
tainty in the value & or to a 40% uncertainty in the median 
redshift of the samples. 

The observational data not included in Table Q is in 
good agreement with the values in the table, but is excluded 
to avoid redundancy. The agreement of the redundant data 
provides further confirmation and indicates that errors are 
dominated by sampling variance rather than by the method- 
ology or the systematics. 



4 RESULTS 

FigEl compares the wtg observations with predictions for a 
fixed value of Q m = 0.3 and three different values of Qa ■ 
We can see how the shape of the prediction depends on the 
amount of dark energy. Even though Wisw at z = depends 
only weakly on Qa , the evolution with redshift depends more 
strongly on Qa- For a fixed Q m , models with larger values of 
Qa evolve more rapidly with redshift to the EdS case, where 
the ISW effect vanishes. Thus, contrary to what happens at 
z = 0, the lower the value of Qa (for a fixed Q m ) the larger 
the ISW amplitude at high redshifts (see also Fig^J- 
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ISW cross-correlation with WMAP 



2MASS 




0.1 1 

mean redshift (z) 
Figure 4. Symbols with error bars correspond to the different 
measurements myg/i) in Table 1. As an illustration of the shape, 
the continuous, short-dashed and long-dashed lines show the con- 
cordance (n m = 0.3,0a = 0.7), opened (f2 m = 0.3, Ha = 0.0) 
and closed (Q m = 0.3, SIa = 1-1) model predictions (at 9 = 6°). 
The dotted line corresponds to the galaxy-galaxy prediction (and 
also the dust contamination model). All lines have arbitrary nor- 
malization. 

To test model predictions with the data, we use a stan- 
dard x 2 -test, x — ^\ (Oi — Ti) 2 /a 2 , where Oi and at 
correspond to the different measurements and errors and T 
correspond to the model. The label i runs for i = 1 to i = 5 
marking the different data points (column 1 in Table as 
we move in redshift. In order to take into account the error 
in the median redshift we take: 

2 ,2 . ( d(w tg /b) \ 2 a 
°i = <*w + I T" I °z (11) 

where a z and a w are the errors in the and z re- 

spectively (see Table 0. We use the relative \ 2 values, 
X 2 ~ Xminy to define confidence levels in parameter esti- 
mation. Top panel of FiglBI shows the resulting confidence 
contours. Taking T = we evaluate the significance of the 
combined ISW detection. We find that this null hypothe- 
sis is rejected with a very high probability: P ~ 99.997% 
(from P„= 4 (x 2 > 26) ~ 3 x 10~ 5 ). We next compute the 
expected ISW effect and compare it with the observational 
data within the ACDM family of models, where O m , Oa 
and h 2 are free parameters (we fix the baryonic content 
fib ~ 0.05 and the primordial spectral index n 3 ~ 1). We 
choose to normalize each model by fixing as- To make our 
results independent of this normalization we will marginal- 
ize as over the range as = 0.8 — 1.0 (flat prior used). As we 
compare u>gt normalized to the CCM model bias, we need 
to compute the relative bias for other LCDM models (see 
section H2.2H . We choose to estimate this at R — 8Mpc/h, 
but the actual number has little effect in the conclusions. 
We have also marginalized over h in a flat prior range 
h = 0.72 — 0.77). Our results are not very sensitive to the 

2 We use H = 100 h km/s/Mpc. 
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Figure 5. One, two and three sigma confidence contours in the 
(fi m , Oa) plane (marginalized over h) for the ACDM model. Top: 
constraints from only ISW. Bottom: constraints from SNIa (blue) 
and ISW (green) along with the combined contours (purple). 



ranges used for as and h: increasing these ranges by a factor 
of two change our contours in less than 20 %. 

The best fit using only ISW data corresponds to £l m ~ 
0.26 ± 0.08, Ov ^ 0.82 ± 0.40, in good agreement with 
other cosmological probes mentioned above. Bottom panel 
of Fig. refcombinedconts we show the confidence contours 
for a ACDM model along with the constraints from recent 
SNIA (Barris et al. 2004) observations. From the figure it 
is clear how the ISW effects gives new complementary infor- 
mation about the cosmological parameters. The EdS model 
is ruled out to high significance. The confidence contours 
are almost perpendicular to the SNIA contours, allowing to 
constrain the parameter space of the model well with just 
these two observations. Combination of ISW with supernova 
data yields fi A = 0.71 ± 0.13 and fi m = 0.29 ± 0.04. 



4.1 Uncertainties in the selection function 

We explore here how robust are our results to the uncer- 
tainties in the galaxy selection function. We take a generic 
parametric form of the type: 

<K*)d* = ^^P^ei-^y 'dz (12) 

so that it is normalized to unity. Parameters /3 and m 
control the shape of the function and are treat as fix param- 
eters; zo is being changed accordingly to the median redshift 
z we want for the selection function. When computing our 
results we use (3 = 1.5 and m = 2, in which case z = \A\zo- 
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— P=l.5m=2 
- |i=2.5m=4 




Figure 6. Two different selection functions with the same me- 
dian redshift z = 1.41. Both have the generic form given by the 
equation 1121 . The black continuous line corresponds to /3 = 1.5 
and m = 2 while the red dashed line is for /3 = 2.5,m = 4 



In order to clarify the role of the selection function 
shape we recalculate our results with a much more peaked 
selection function. This second selection function have (3 = 
2.5 and m = 4 and it is plotted together with the fidu- 
cial one in Figure H Both cases have the same median red- 
shift z — 1.41. Top panel of FiglTlshows the contours in the 
($l m , Qa) plane for the more peaked selection function (with 
P = 2.5 and m = 4). The contours are similar to the fidu- 
cial model (ie compare to Fig. |SJ but favoring slightly lower 
values for Qa and fl m . 

Besides the uncertainty on the shape of the selection 
function there is also uncertainty in the median redshift. 
We have checked what happens if this uncertainty is not 
taken into account. We just set the redshift errors a z = 
in Eg. 1111 . Contours for the (fi m ,fiA) plane are plot in the 
bottom panel of figure[7| which are also to compare with fig- 
ure |S] There is hardly any difference because the theoretical 
values of wtg change very little within the median redshift 
error range. 



4.2 Equation of state 

The ISW effect can also be used to constrain the dark en- 
ergy equation of state parameter. In this case, as suggested 
by the CCM, we assumed a flat universe. We focus on a 
constant w parameter and maintain the same flat priors for 
h and <t 8 ( 0.72 < h < 0.77 0.8 < a 8 < 1.0 ). Top panel of 
Figure |H| show the one,two and three sigma contours for the 
(Q m ,w) plane using only the ISW data. Join contours with 
the SNIa data are shown in the bottom panel of Figure |S] 
Both datasets are also complementary for the w determina- 
tion. The SNIa data is from (Barris et al. 2004). 

Making a join ISW+SNIa analysis with the flat prior 
reduces notably the allowed space for the parameters to 
w = -1.02 ± 0.17 and Q A = 0.70 ± 0.05. The contours are 
comparable with other analysis in literature (Sandvik et al 
2004) which combines SNIa with WMAP and SDSS data. 




0.25 0.5 




0.25 0.5 



Figure 7. One, two and three sigma confidence contours in the 
(H m ,OA) plane for the ACDM model. Top panel: contours using 
a more peaked selection function (/3 = 2.5, m = 4) but with the 
same mean redshift as the fiducial case (ie compare to Fis;l5l. 
Bottom panel: contours when errors in the median redshift of the 
selection functions are neglected. 

The results we found are still in full agreement to the CCM 
with S7a ~ 0.7 and w = — 1. 



4.3 Possible Contaminants 

The constraining power of the new ISW data comes from 
the simultaneous fitting of data at different redshifts, that 
is from the shape information in Fig. [I] Because of the uncer- 
tainties in the relative normalization due to a relative bias, 
any given point alone does not constrain well the cosmolog- 
ical parameters. But the combination of the data gives us a 
new powerful tool for cosmological parameter estimation. 

The shape of the curve as a function of redshift also 
provides an important test for systematics. CMB and galaxy 
maps are both masked and corrected from galactic absorp- 
tion/extinction, but any residual contamination could pro- 
duce a cross-correlation signal. Emission and absorption by 
our own galaxy produce patchy hot spots in the CMB maps 
and negative density fluctuations in the galaxy distribution 
(because of extinction). In principle, this should therefore 
result in a negative cross-correlation, but overcorrecting for 
the effects of galactic absorption could also result in a pos- 
itive signal. This possibility have been tested for each of 
the samples, by comparing the cross-correlation to WMAP 
maps at different frequencies. Most analysis use the WMAP 
KpO mask, which excludes about 30% of sky on the basis of 
galactic or extra-galactic (eg radio sources) contamination. 
In all cases the contamination seems smaller than the errors 
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Figure 8. One, two and three sigma confidence contours in the 
(fiA,if) plane (marginalized over h and ag. Top: constraints only 
from ISWfrom ISW. Bottom: constraints from SNIa (blue) and 
ISW (green) along with the combination (purple). 

(eg see Fig. 2 in Fosalba & Gaztanaga (2004)). Moreover, 
one does not expect this effect to have any redshift depen- 
dence, contrary to the measurements in Fig. [I] 

Cold dust in distant galaxies, will also produce patchy 
hot spots in the CMB maps and positive density fluctuations 
in the galaxy distribution (could also be negative because of 
internal extinction). The resulting cross-correlation should 
trace the galaxy-galaxy auto correlation function, wgg, and 
should therefore have a very different redshift dependence 
to the ISW effect. The dotted line in Fig. 0] shows the pre- 
dicted shape dependence for wgg contamination with arbi- 
trary normalization. The shape is clearly incompatible with 
the actual cross-correlation measurements. It is also worth 
noting how wgg goes quickly to zero at z ~ 0.2, while the 
ISW cross-correlation remains positive. This is due to the 
fact that at these corresponding large scales, > 40 Mpc/h, 
matter-matter correlations w mm effectively decays to zero, 
while wtg, which traces the gravitational potential, has a 
less rapid decay with distance. 



5 CONCLUSION 

The cross-correlation of CMB anisotropics with very dif- 
ferent galaxy surveys provides consistent detections. Their 
combination follows the CCM predictions with a probability 
of only ~ 3 x 10~ 5 for being a false detection. This provides 
new and independent evidence for dark energy and dark 
matter, ruling out the EdS model to a high significance (for 
any value of Ho)- Combination with SNIA data results in 



strong constraints to ft A = 0.71±0.13 and Q. m = 0.29±0.04. 
This in good agreement with the flat universe Q m + f^A — 1 
found independently by CMB data (Bennett et al. 2003; 
Tegmark et al. 2004). If we assume a flat universe, we find 
Qa = 0.70 ± 0.05 and w = -1.02 ± 0.17 for a constant dark 
energy equation of state. The data shows, for the first time, 
statistical evidence of a recent slow down in the growth of 
structure formation on linear scales, just as expected in a 
flat accelerated universe. The new ISW constraints rely in 
a totally different physical effect that previous cosmological 
constraints, providing new light on a dark cosmos. 

Note added in proof: After this paper was originally 
submitted to astro-ph ( astro-ph/0407022 ) a related analysis 
using our data compilation have been published by Corasari- 
niti, Giannantonio and Melchiorri (astro-ph/0504115 ). 
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